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Goal  of  this  project 

Investigations  of  arbitrary  complex  functions  for  electromagnetic  waves  by 
using  a  plasma  array  reinforced  by  metallic  components  and  changing  a  pattern  of  series 
switches,  external  power  for  plasma  generation,  and  pressure  of  discharge  gases. 

I.  Introduction 

A  metamaterial,  which  has  a  functional  microstructure  leading  to  an 
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extraordinary  macroscopic  property,  is  composed  of  metals  and  dielectrics  in  general.  ' 
Specific  and  careful  design  of  its  spatial  structure  is  required  before  it  works,  for 
instance,  as  a  negative-refractive-index  (negative-  N  )  material. 4  (1  To  secure  flexibility 
of  its  design  and  parameters,  here  we  propose  plasmas  as  materials  composing 
metamaterials.  Microplasmas,  which  have  sizes  ranging  from  millimeter  to  micrometer 
scale  and  work  near  and  at  atmospheric  pressure,  have  electron  densities 
ne  =  1012 -TO16  cm  3  that  correspond  to  electron  plasma  frequencies  copsl 2n  at  10 
GHz-1  THz.  They  are  equivalent  dielectrics  with  permittivity  s  between  unity  and 
zero  above  cope ,  and  equivalent  metals  with  negative  s  below  a>pe .  ne  and  the  gas 
pressure  that  determine  s  are  controllable  by  adjusting  input  of  electric  power  and  gas 
inlet  for  plasma  generation,  and  their  spatial  profiles  can  be  altered  to  a  certain  extent. 

We  have  reported  plasma  photonic  crystals  that  exhibit  significant  functions  of 


an  array  of  microplasmas  for  control  of  electromagnetic  waves.7'10  In  a  plasma  photonic 
crystal,  a  1 -dimensional  or  2-dimensional  array  of  spatially  periodic  plasmas  is  fonned, 
and  topical  features  such  as  photonic  band  gaps  and  surface-plasmon-like  chains  were 
predicted  theoretically  and  observed  experimentally.  We  also  pointed  out  the  importance 
of  complex  a ,  which  can  be  controlled  in  a  discharge  plasma  and  makes  a  band  gap 
around  which  wave  attenuation  changes  drastically.  In  such  a  structure,  £  was 
elaborately  controlled  in  a  homogeneous  array,  although  permeability  //  was  constant 
at  unity. 
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Here,  we  extend  these  results  to  metamaterials  composed  of  plasmas, 
referred  to  as  "plasma  metamaterials"  hereafter.  The  concept  of  plasma  metamaterials 
was  proposed  in  Ref.  1 1 .  One  type  of  plasma  metamaterial  was  shown  in  Ref.  12,  which 
verified  that  spoof  surface  plasmons  propagating  on  a  perforated  metal  plate  can  be 
controlled  dynamically  by  microplasma  generation  in  the  holes  of  the  plate.  Another 
type,  possibly  with  the  most  significant  role  that  plasma  metamaterials  may  play,  is  a 
dynamic  negative-refractive-index  material  composed  of  negative-  e  plasmas  and 
negative-  //  (metallic)  structure;  our  brief  previous  report  demonstrated  experimental 
results  of  plasma  metamaterials  with  negative  refractive  index,  although  theoretical 
background  has  not  been  clarified  thus  far.  Theoretical  design  of  plasma  metamaterials 
is  required  for  further  extension  of  plasma  metamaterials  in  industrial  applications  as 
well  as  for  scientific  purposes. 

Before  we  describe  plasma  metamaterials  specifically,  we  take  a  close  look  at 
their  performance  requirements.  When  we  generate  plasmas,  metallic  discharge 
electrodes  are  usually  required.  If  we  design  metallic  wires  as  electrodes  to  generate 
microplasmas  with  negative  a  and  simultaneously  as  resonators  to  make  macroscopic 
//  negative  due  to  magnetic  resonance,  the  synthesized  structure  with  spatial 
periodicity  is  expected  to  be  a  state  with  negative  A  as  a  double-negative  system.  We 
note  that  the  state  with  negative  a  arises  from  microscopic  parameters  in  plasmas, 
which  differs  substantially  from  macroscopic  permittivity  emerging  in  an  array  of  thin 
metal  wires.14'15 

Another  important  characteristic  of  plasma  metamaterials  is  the  role  of  the 
imaginary  part  of  a , 10  at ,  given  as 
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This  is  mainly  associated  with  the  electron  elastic  collision  frequency  with  neutral 
particles  vm ,  which  directly  depends  on  gas  pressure.  ax  is  estimated  as  a  loss  source 
for  propagating  electromagnetic  waves,  and  we  usually  tend  to  minimize  its  value  to 


avoid  wave  attenuation.  Figure  1  shows  complex  values  of  (sr,  in  various  gas 
species  and  pressures  with  several  levels  of  nt ,  which  indicates  that  we  can  control  ( st , 
sx )  freely  by  changing  gas  conditions  and  energy  input  for  generation. 

However,  when  we  consider  the  refractive  index  N  =  Nr  +  j Nt  as  a  product  of  s  and 
M,  say,16 

N  =  y[s^u,  (2) 

we  recognize  that  wave  attenuation  expressed  by  Nx  is  not  directly  connected  to  sx . 
To  understand  wave  propagation  more  rigorously,  calculation  of  a  dispersion  relation  in 
the  periodic  structure  with  complex  N  is  required,  where,  in  the  dispersion  relation, 
wavenumber  k  =  k1+jki  is  as  a  function  of  electromagnetic  wave  frequency  col  2n. 
In  our  previous  report  where  we  mentioned  a  specific  plasma  metamaterial  with 
non-one  //  using  experimental  results,  we  assumed  that  both  a  and  //  are  real 
numbers  for  simplicity,  which  is  not  sufficient  for  understanding  the  importance  of 
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In  this  project,  we  investigated  arbitrary  complex  functions  for  electromagnetic 
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Fig.  1.  Complex  e  as  functions  of  gas 
species,  pressure  and  ne  deduced  from  Eq. 
(!)• 


II.  Theoretical  prediction  of  flexible  plasma  metamaterials 

Detailed  formulation  of  the  theoretical  model  was  described  in  Ref.  17.  In  short, 
Faraday’s  law  and  Ampere’s  law  are  applied  to  calculate  an  equivalent  // .  Relationship 
between  induced  voltage  and  current  is  deduced  from  Foster’s  theorem.  And  effects  of 
micro-periodicity  whose  length  is  much  less  than  a  wavelength  is  included  with  the 
knowledge  of  the  Bloch  modes. 

In  the  following  results,  we  used  Eq.  (2)  for  derivation  of  N  where  s(=  e  ) 
is  expressed  by  Eq.  (1)  and  displayed  in  Fig.  1,  and  //(=  which  is  a  pure  real 
number,  is  displayed  in  Fig.  2.  k  is  derived  by  solving  Eq.  (8),  where  kr,  the  real  part 
of  k ,  has  the  maximum  value  of  12.57,  which  corresponds  to  the  end  of  the  first 
Brillouin  zone. 

Figure  3  shows  frequency  dependence  of  N  and  dispersion  relation  when 
collisionless  plasmas  with  apx  =  0  are  generated.  Here  we  set  ne  =4x10  "cm'  to 
derive  sp ;  in  the  following  calculation,  we  use  this  ne  value.  spx  is  always  negative 
at  frequencies  of  less  than  1 8  GHz  where  co  =  a>pe ,  and  so  N  becomes  imaginary  at 
most  frequencies  outside  the  resonance-vicinity  bands  and  wave  propagation  is 
prohibited.  However,  since  both  sp  and  //eff  are  negative  in  the  resonance-vicinity 
bands,  N  becomes  real  and  negative.  As  a  result,  kr  has  a  finite  value  while  k]  is 


Fig.  2.  Macroscopic  m  for  double-helix 
metal  wires  with  calculated  admittance  Y 
for  structure  in  Fig.  2.  ID  plasma 


Fig.  3.  Frequency  dependence  of  N  and  k  of 
ID  collisonless  plasma  metamaterial  with 
ne  =  4xl012  cm'3. 


zero  in  these  bands,  and  waves  can  propagate  with  negative  phase  velocity.  Arising  from 
the  fact  that  N  <  0  ,  several  extraordinary  phenomena  are  expected,  such  as 
performance  as  a  perfect  lens." 

Now,  we  proceed  further  to  look  at  the  effects  of  finite  spi .  Figure  4  shows 
frequency  dependence  of  N  and  dispersion  relation  when  plasmas  are  slightly 
collisional  at  50  Torr  of  Ne.  In  the  resonance-vicinity  bands,  Nr  is  negative,  similar  to 
the  value  in  Fig.  3,  but  Nt  is  not  zero.  Consequently,  kt  is  a  certain  value  while  kr 
is  still  larger;  the  waves  propagate  with  slight  attenuation.  We  note  that  spx  is  so 
negative  that  we  easily  recognize  that  N r  is  negative  from  Equation  (1).  Figure  5 
shows  frequency  dependence  of  N  and  dispersion  relation  when  plasmas  are  more 
collisional,  at  150  Torr  of  Ne.  In  the  resonance-vicinity  bands,  Nt  is  still  negative,  but 
Ni  is  large  with  its  amplitude  comparable  to  Nr .  Consequently,  k]  is  almost  similar 
to  kr .  In  this  case,  spr  is  negative  with  larger  spi  value. 

Figure  6  shows  frequency  dependence  of  N  and  dispersion  relation  when 
plasmas  are  highly  collisional  at  760  Torr  of  He.  In  the  resonance-vicinity  bands,  Nr  is 
still  negative  as  well,  but  now  Af  is  much  larger  than  that  in  Fig.  5.  Consequently,  C 
is  quite  large,  although  kr  is  still  finite.  We  have  to  point  out  that,  in  this  case,  spr  is 
not  negative,  but  Nr  is  certainly  negative.  This  phenomenon  is  unique  to  a  plasma 
metamaterial,  in  which  spi  plays  an  important  role.  When  we  take  a  careful  look  at  Eq. 
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Fig.  4.  Frequency  dependence  of  N  and  k  of 
ID  lossy  plasma  metamaterial  with  ne  = 
4x10 12  cm'3  at  50  Torr  of  Ne. 


Fig.  5.  Frequency  dependence  of  N  and  k  of 
ID  lossy  plasma  metamaterial  with  ne  = 
4xl012  cm'3  at  150  Torr  of  Ne. 


(1)  again  with  possibilities  of  complex  a  and  //,  we  recognize  a  possibility  that  large 
api  with  respect  to  small  apr  may  contribute  to  yielding  negative  Nr  with  large  Nt . 

From  these  results,  we  can  recognize  that  k]  and  N{  can  be  set  to  different 
kr  and  Nr  values  by  changing  gas  conditions.  As  we  mentioned  earlier,  Nr  and  Af 
affect  phase  shift  and  amplitude  attenuation  of  waves,  and  so  we  can  predict  potentials 
for  complex  value  filters  by  plasma  metamaterials. 

III.  Experimental  verification  of  flexible  plasma  metamaterials 

1  'X 

We  briefly  reported  a  controllable-  N  material  in  preliminary  experiments. 
Double-helix  metal  wires  and  generated  plasmas  formed  this  metamaterial.  In  the 
double-helix  structure,  inductance  L  is  expected  to  be  along  a  wire,  and  capacitance 
C  is  present  between  the  wires.  L  and  C  form  a  ladder-like  equivalent  circuit, 
showing  a  series  resonance  when  two  of  its  ends  are  connected.  Consequently,  the  entire 
structure  has  multi  resonance  frequencies  above  which  macroscopic  ju  of  its  array 
becomes  negative.  As  shown  in  Fig.  7,  the  array  of  the  double-helix  metal  wires  was  set 
in  a  capillary  tube,  and  the  tubes  were  installed  on  a  coplanar  microwave  waveguide  to 
investigate  the  entire  structure  as  a  device  under  test  (DUT).  Interference  measurement 
indicated  that  the  DUT  became  negative-  N  material,  although  the  results  obtained  in 
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this  measurement  include  large  error  bars. 

To  confirm  negative  N  in  this  metamaterial  configuration  in  a  more  precise 


Fig.  6.  Frequency  dependence  of  N  and  k  of 
ID  lossy  plasma  metamaterial  with  ne  = 
4xl012  cm'3  at  760  Torr  of  Fie. 


Fig.  7.  Schematic  view  of  device  under 
test  composed  of  double-helix  metal 
wires  and  generated  micronlasmas. 


diagnostics,  we  developed  a  microwave  lumped-element  circuit  for  a  parameter  retrieval 
method  using  the  two-port  twelve-term  error  model,  which  in  displayed  in  Fig.  8. 
In  short,  signals  are  analyzed  by  splitting  into  in-phase  and  quad-phase  parts,  and 
precise  phase  information  is  deduced  in  the  assumption  of  sinusoidal  wave  forms,  where 
quad-phase  signals  have  phase  difference  of  90  degrees  from  the  initial  one.  The 
measured  microwave  properties  through  the  DUT  are  partly  shown  in  Figs.  9  and  10, 
where  various  gas  conditions  were  set  to  see  the  effects  of  complex  s  shown  in  Fig.  1. 
First,  the  time  evolutions  of  S21  were  very  different  from  those  of  transmitted  signals 
in  the  2D  plasma  photonic  crystal.  ~  The  composites  of  plasmas  and  metallic  resonators 
make  the  microwave  transmitted  signals  vary  more  drastically  in  both  amplitude  and 
phase. 

Here  we  estimate  values  of  N  in  a  typical  case  shown  in  Fig.  9  where 
metamaterial  effects  are  dominant.  Without  plasmas,  the  resonance  frequency  of  the 
array  of  the  double-helix  wires  was  around  4.25  GHz,  where  Re(N)~1.8  with  the 
length  of  the  structure  =  3.8  mm  and  the  corresponding  phase  shift  through  the  structure 
-  6.1  rad.  At  4.27  GHz,  which  was  very  close  to  the  resonance  frequency,  and  for  Ne  at 
20  kPa,  the  change  of  arg(S21)  by  plasma  generation  was  more  than  2n  rad,  which 

indicates  that  Re  (A/)  of  the  DUT  decreased  down  to  negative  values.  From  calculation 
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of  the  filling  factor  of  the  double-helix  wires  (-0.32)  in  the  entire  DUT  structure, 


Reflected  signal 


Signal 

Generator 


DUT 


Directional  Coupler 
Sig.LQ 


Power  Splitter(A) 


(b) 


Power  Splitter  (B) 


Sig.LQ 


DBM1 


... 

,  F1 

■ 

RF1 


Hybrid  divider 


Sig.  LO 


Oscilloscope 


LOl 


L02 


Power  Divider 


Sig.  RF 


DBM2 


RF2 


(b) 


~o 

<v 

■*->  , 
■*-> 

E 

( A 

Sig.RF 

1— 

Sig.RF 

Fig.  8.  Detector  circuit  for 
parameter  retrieval  method  of 
dynamic  metamaterials. 


Re(N)  in  the  metamaterial  layer  composed  of  the  microplasmas  and  the  double-helix 
metal  wires  was  estimated  to  be  -1.5  -  -3.0. 

If  the  discharge  gas  conditions  are  different,  corresponding  transmitted 
microwave  signals  are  along  different  contours  in  the  complex  S21  plane,  as  shown  in 
Fig.  11.  For  instance,  in  the  case  of  Ne  at  20  kPa  in  Fig.  9,  the  phasor  of  S21  rotated  at 
more  than  2n  rad  since  complex  s  changes  in  both  real  and  imaginary  components,  as 
shown  in  Fig.  1 .  We  also  note  that  this  gas  condition  roughly  corresponds  to  the  case  in 
Fig.  5.  On  the  other  hand,  in  the  case  of  He  at  100  kPa  in  Fig.  10,  the  phasor  of  S21 
swung  along  one  straight  line  which  implies  no  substantial  change  of  the  phase 
component.  The  corresponding  s  might  change  only  in  its  imaginary  component  with 
almost  no  change  of  Re(^),  which  roughly  corresponds  to  the  case  in  Fig.  6.  These 
facts  verify  that  the  complex  components  in  s  and  N  are  key  factors  and  distinguish 
the  plasma  metamaterials  from  the  ordinary  solid  metamaterials. 


IV.  Summary  and  future  perspectives 

We  theoretically  and  experimentally  verify  complex  electromagnetic -wave 
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filters  achieved  by  plasma  metamaterials.  In  addition  to  the  previous  results  in  which 
number  of  turn-on  plasma  columns  was  controlled,  we  can  control  electromagnetic 
waves  with  adjustable  phase  shift  and  attenuation  by  changing  gas  pressures  and 
external  power  supply  for  plasma  generation.  Although  such  a  concept  and  verified 
technology  are  applicable  for  alternatives  of  conventional  solid-state  microwave  devices, 
we  cannot  recognize  certain  advantages  to  replace  them  by  novel  devices  based  on 
plasma  metamaterials.  Instead,  when  we  apply  plasma  metamaterials  to  controllable 
media  for  high-power  microwaves  working  as  electric  energy  carriers,  they  will  be 
superior  to  other  conventional  devices  since  plasmas  do  not  deteriorate  by  microwave 
with  high  energy  density. 


Re(S21) 


Fig.  11.  Time  trajectories  of  complex 
phasor  S2i  in  two  corresponding  cases  in 
Fig.  9  as  (a)  and  Fig.  10  as  (b). 
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